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ABSTRACT 

The true and eccentric anomaly parametrizations of the Kepler motion are 
generalized to quasiperiodic orbits, by considering perturbations of the radial 
part of the kinetic energy in a form of a series of negative powers of the orbital 
radius. A toolchest of methods for averaging observables as functions of the 
energy E and angular momentum L is developed. A broad range of systems 
governed by the generic Brumberg force and recent applications in the theory 
of gravitational radiation involve integrals of these functions over a period of 
motion. These integrals are evaluated by using the residue theorem. In the 
course of this work two important questions emerge: (1) When does the true 
and eccentric anomaly parameter exist? (2) Under what circumstances and why 
are the poles in the origin? The purpose of this paper is to find the answer to 
these queries. 

Subject headings: binaries: close, gravitation, pulsars: general — relativity 
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1. Introduction 



In a recent series of papers ( [Gergely et al., 1998a[ - |Cergely et al., 1998c|) , the present 



authors have derived the cumulative back reaction effects of gravitational radiation emitted 
by a binary system of spinning components. The purpose of our work has been to study the 
effects of the spins on the secular evolution of the system under gravitational back reaction, 
necessary in order to provide realistic signal templates for the designated gravitational 
wave observatories LIGO ( [Abramovici et al., 199!^ ), VIRGO ( [Bradaschia et al., 1990D , 



GEO ( [Hough 1992| ) and TAMA ( |Kuroda et al., 1997| ). Neutron star-black hole and black 



hole-black hole binaries are designated sources (Thorne 1997) for the Laser Interferometer 



Space Antenna (LISA, Pender et al., 1996| , [Schutz 1997|) . In the course of the computation. 



two parameters of the quasiperiodic orbit have been efficiently utilized for evaluating 
integrals of the type 

rT ,.j 

-dt , (1^ 



where T is a radial period and u; is a function that has no manifest radial dependence (in 
a sense that will be made more precise in Section III). The radial period, for instance, is 
determined by this integral with n = —2, while the instantaneous radiative losses of the 
quantities characterizing the orbit are sums of contributions of type (|l]) with n = 1, .., 7. 

For averages with > 0, a parameter x generalizing the true anomaly of the Kepler 
motion renders the poles of the integrand in the origin (after passing to the variable 
z = e*^), thereby easing the evaluation of the integral. For n < 0, on the other hand, the 
parameter ^ generalizing the eccentric anomaly proves suitable. Hence in both cases the 
residue theorem can be applied to integrands with a pole in the origin of the complex plane 
and, in certain cases, with an additional pole. This is to be contrasted with the elaborate 
structure of the integrands when performing a radial integration over the same quantities 
( IRleth t Schafer l997| ). 
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Our purpose in this work is to gain an insight into the nature of these two 
parametrizations of the perturbed Kepler motion, and to systematically search the structure 
of the integrals of the type (|l|). We show that these parametrizations exist for a wide class 
of quasiperiodic orbits. Our method of averaging is independent of the problem of radiation 
reaction. We therefore expect that the method developed here will find applications, beyond 
the theory of gravitational radiation, in the treatment of astrophysical binary systems. 

We will study perturbations of the Kepler motion represented by a series of negative 
powers of the radial coordinate r in the expression of r^. This type of perturbation is 
generic enough to include both the radiation reaction perturbative terms of our previous 
works as well as to comprise other physically relevant situations. 

A wide class of perturbed Kepler motions, which fit in our description, is described by 
the Brumberg force ( Brumberg 1991 , [Sottel 19"S5 ) 



mr mr 



'2'ym ^ 9 (rr)^l 2Am(rf) . 



(2) 



which can be derived from the Lagrangian 



(3) 

Here a, /5, 7 and A are perturbation parameters, the total mass is m = mi + m2 and 
the reduced mass is /x = 17111112/ m. 

The energy E = (dC/dr) r—C and the magnitude L of the angular momentum 

^ ^ / 2a -2/3 + A ^ 2Am\ 

L = Ln 1 + E + (4) 

V /i r J 

of the system are conserved. In (|^) Ln = yur x r is the Newtonian angular momentum. By 

substituting in the Keplerian part of the energy E the identity 

= + 4^ (5) 
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(with Lj^ expressed from (Q) in terms of L^) and inserting the Keplerian expression 
= 2E/fi + 2m/ r in the perturbative terms of E, the radial equation of motion is found: 

r2 = — + — - {2a~2p + \)—-~{3a-2p + 2\) (6) 

jj, r /i^r^ /i"^ /ir 

2ELF' 2TnIj^ 
- (3a - 2/5 + 2A + + (2a - 2/5 + A)^^ + (a + 2A)^-^ . 

With a = Z//2, /? = (! + 3i/)/2, 7 = 2 + i/, A = 2 — i/we get the radial equation of 
Damour fc Deruelle, 1985| . 



The other relevant example which fits in the generic perturbation scheme to be treated 
in this work is the dynamics of the spinning binary system ( [Barker fc O'Connell, 198^ , 



Kidder 1995| ). The major difference compared to the Brumberg-type perturbation is that 



the orbital angular momentum L is not parallel to Ln. For comparison, we give the radial 
equation characterizing the spinning binary system, up to linear terms in spin (|Gergely~et 
|al., 1998c|) : 

, 2E 2m 2E(L-a-) 2(2L ■ S + L-ct) 
r = 1 \ ^ . (7) 

Here S and cr are the total and the weighted spin vectors: 

S = Si + S2 , cr = r/Si + r7""^S2 , 1] = — . (8) 

mi 

In the next section we study the radial component of a perturbed Keplerian motion 
generic enough to include both (||) and (|^). We express the turning points of such a radial 
motion in terms of constants of the motion. 

Then in Sec. Ill, we solve the radial motion by introducing the true anomaly 
parametrization and we formulate rigurous statements about the poles of the integrands 
(|I|), for n > 0. Integrands with n < will be discussed in a similar fashion in Sec. IV, 
after introducing the eccentric anomaly parametrization. An immediate application of the 
developed methods is the computation of the radial period for a wide class of perturbed 
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Keplerian motions. This is done in Sec. V. Finally in Sec. VI. we summarize our results. 
Throughout the paper we use relativistic units G = c = 1. 



2. The radial motion 



We consider the radial equation of a perturbed Kepler problem|^ 

where the small coefficients couple the perturbations to the system. The ffist three terms 
on the right are familiar from the radial Kepler motion, however, the constants E and L 
characterize the perturbed motion. In a Lagrangean treatment, the conservation of the 
energy E and the total angular momentum J follows from symmetry considerations. We 
assume these symmetries to characterize the system to the required order. Furthermore, 
we assume that the total spin S of the system evolves in a direction orthogonal to the 
orbital momentum L, a requirement which is met by the spinning binary system and has 
the immediate consequence that L is conserved. In what follows, we shall employ a zero 
subscript for denoting the parts of expressions taken from the Kepler motion. For example, 

2 2m 

^0 = — + ^ ■ 10 

The turning points at = are the solutions of the algebraic equation: 

p 

2^ErP + 2m/iVP-i - LV^'^ + ^ ip^rP-' = 0. (11) 

i=0 

Without the perturbing terms this Eq. has the roots: r+, r_, 0,..., 0, where there are 



^ For a review of the Newtonian and relativistic perturbed Kepler problem, cf. p amour 



fc Deruelle, 1985 
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p — 2 zero roots 



-2E /i(m/i^Ao) ^ ^ 



and Aq is the length of the Runge-Lenz vector, given by 



(13) 



Then the roots of Eq. (|T^) describing the perturbed motion can be found in the form: 
Tinax = r+ + rinin = r_ + £_ and (for p > 2) e^, k = 1, ...,p — 2, where Ek is real or 
imaginary: 



£k 



l/{p-2) 



(14) 



and 



£± = - 



+ 2m/i^(p - l)r± - (p - 2)1.2] ^ ^.^3- 

j=0 



3-i 



(15) 



The roots e^, even if real, are unphysical as they are smaller than the radius r. The 
exponent l/{p ~2) renders this approximation increasingly corrupt as p ^ 00. 



Thus the turning points are: 



1 ^ 



min —2E 2fiAo -^Q 

We will need also their reciprocal: 



/i(m/i =F ^0) 
L2 



i-2 



(16) 



1 /i(m/i =F Aq) 1 



max 
min 



L2 



fi{mfi =1= 

Z2 



2/iAo ^ 

Note that is defined by Eq. (pTsP, with the constants and L characterizing the 
perturbed motion. 



(17) 
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3. True anomaly parametrization 

The true anomaly parameter x of the Kepler ellipse r = p (1 + ecosx)"^ lying in the z 
plane with semilatus rectum p = and eccentricity e = — is the azimuthal angle. A like 
parametrization can be introduced for the perturbed motion, but this can be done in several 
inequivalent ways. An example of this parametrization has been discussed by P^yan 19% 
for the radial motion of a test particle orbiting a spinning point mass. When one attempts 
to utilize his parameter for averaging the radiative losses, the integrals involved will contain 
poles at various location in the complex x plane; and this feature is not amenable for doing 
the integration. However, adding any term proportional to sin x will not alter the turning 
points in r. One may take advantage of this freedom in fixing the generalized true anomaly 
parameter in such a way that two main features, one relating the turning points with the 
scale of the parameter and the other the simple periodic dependence of the radius on the 
parameter, are stipulated for the perturbed motion: 



(a) rrain = r{0) 

rmax = r{n) (18) 

dv 1 

(b) — = -7r2 - = 7COSX + 5 ■ 

a(cosxj r 

These two properties are equivalent with: 

2 ^ 1 + COSX ^ 1 - cosx .^g^ 

' ' mm ' max 

The true anomaly parametrization has a remarkable feature, which can be encapsulated in 

Theorem 1: The integral J ■:pi+^dt , n > is given by the residue in the origin of the 
complex true anomaly parameter plane z = e^^. 



- 9 - 



In order to establish this theorem, first we pass in the integral to the true anomaly 
parameter: 



1 



dt 



1 dt dr ^ 1/1 1 
-dx 



mm ' max 



1 sinx 



dx- 



(20) 



j.2+n 2 Vr„ 

Here we have used (plsl) . The last integrand is a product of two factors. The parametric 
form of the first one, is given by the binomial expansion 

,j /I + cosx\* /I — cosx' 



n n 



i=0 



(21) 



Hence this factor is a polynomial in cos x- 



Next we will show that the second factor sinx/r is also regular. For this purpose we 
write the radial equation of motion (^ ) in terms of the parameter x- Using (|2l|) , the Kepler 
part (|10D of becomes: 

" ^ 1 1 2mu^ / 1 1 



A2 t2 
■2 _ ^0 ^ 

° L2 4/x2 



mm ' max 



2mu^ / 1 



cosx 



max ^ \ ' mm ' m 

We insert here the expressions ( p^ of the turning points r^i^ and rmax in terms of the 
constants of the motion and introduce the notation 



(22) 



Ai = (m/i + AoY ± {mfi - Aq) 



(23) 



Keeping terms in Tq to linear order in (pi, we get 



• 2 ^0 • 2 

'"o = sm X- 



1 A //i^'- 



L2 



A*l + A^cosx 



cosx • 



2/^' .=0 

Adding the first-order terms present in the radial Eq. (P) we finally obtain 

^2 



(24) 



ism^X-gp(^)K.(x), 



(25) 



where we denote 



A*_ cos X + A'^(_ cos^ X 



- {mfi + Aq cosx)* 



(26) 
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A power series expansion in the small coefficients leads to 

+ ^773—2- 1^-2 (ti) Mx)- (27) 



r Ao 2Ai]sin2x^/xnL2 
We will show, that the numerator of the second term is always proportional to sin^ x- 
Indeed 



^.(x) = ECf(m/ir%'S,(x) , (28) 

k=0 

where we denoted 

1 — (—1)'^ 1 + (—1)*^ 
2fc(x) = ^-^cosx + cos^ X - cos'^ X- (29) 

For k = 0, a. sin^ x factor can be taken out. We then find that this property continues to 
hold for all even and odd values: 

So = - sin^ X, Si = S2 = , 
S2(,+i) = cos2x-cos^+''x = sin2xEcos2('^+^)x, (30) 

n=0 

e-1 

= COS X - cos^+^^ X = sin^ X J2 cos^"+^ x • (31) 

n=0 



Therefore no x dependence remains in the denominator of the integrand (27). We have 
thus proved that sinx/r is regular. 

When passing to the complex parameter 

z = e'^, (32) 

the denominator of the integrand (a polynomial in sinx and cosx) will be an integer power 
of z, leaving the origin as the only location of the pole. 

Theorem 1 being proven, we can now specify which are the functions in (|T]), which 
will not destroy the advantageous property of the integrand. As is transparent from the 
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above proof, such a function can he any polynomial of sin x cos x- In particular it can 
have the form 

u = uo + (p{ip,r) , (33) 

where uq is some constant, and the perturbation term ip{ip, r) is a polynomial in sin-^, cos?/' 
and r. As the azimuthal angle ip and the rate of change r of the radius are needed here to 
Keplerian order: 

V' = ^0 + X , = ^ sin X, (34) 

the function is a polynomial in sin x and cos x- 

The results of this section can be summarized as follows. For the perturbed Kepler 
motions (|^), any integral (|I]) with > satisfying (P3[), by use of (|T^, ( P^D and ( pi[ ) can 
be turned into a complex integral in the variable z. This integral is simple to evaluate with 
the residue IZ lying in the origin: 



— at = 2m 7i[z = 0) . (35) 



4. The eccentric anomaly parametrization 

Another convenient parametrization of the Kepler ellipse r = a(l — ecos^) with 
semimajor axis a = = is given by the eccentric anomaly parameter C,- In terms of 
the turning points, the eccentric anomaly parametrization may be defined as follows, 

2r = (1 + COsO'^min + (1 - COsO'^max ■ (36) 

Hence 

dv 

2^ = (rmax - ^min) siu ^ . (37) 
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For a perturbed Kepler motion (||), we use ( p6D as the definition of tlie generalized 
eccentric anomaly parametrization. This parametrization has a remarkable property the 
main consequence of which is contained in the following 

Theorem 2: The integral J ■:^:^dt with n < is given by the sum of the residues in 
the origin of the complex parameter plane w = e*^ and [for perturbative terms with negative 
powers i > 2 — 72 of r in the radial equation (P)] at 

We prove this theorem by passing to the true anomaly parameter in the integral. For 
convenience we introduce n' = — 1 — n>0. 



1 dt = [r-'~^^±d^ = ^JIlELllIlB^ [r^'^^dC (39) 



j,2+n J dr 2 J rr 

In the last equality we have used (pT]). We insert the binomial expansion of (|36D in 
the first factor of the last integrand. In what follows, we will show that the second factor 
sin ^/rf is regular. 



Multiplying the radial equation (g) by r we get 

\2 



rr 



—r^ + 2mr - — + V . 40 

Hence, using again Eq. (^), we get the equation of the associated Kepler motion in terms 
of the parameter ^: 

E 2 

('"'")o 2^ [('"mill '"max) COS -|- '"min ~l~ '"max] ~l~ ^(^^min '"max) COS -|- "^(^min ~l~ '"max) ■ 

(41) 

By simple algebra, we can write this as 



^rr)l = + — 

^ -2fxE 2/i 



'"min ~l~ '"max ~l~ ^ ~l~ ('"min '"max / COS C, 



2 



(42) 
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By use of the quantities introduced in (p3D, and following the steps in the derivation of 
Eq. (p5D , the full radial equation ( ^0| ) is written as 



[rr] 



-2/iE 



p 



Sin e - L t:^ 



i=0 



L'^imjji — Aq cos^) 



j-2 



= (m/i-AoCOS^ 



^ Ai cos ^ + A^+ cos^ ^ 



-2EL'^V 



(43) 



2 \ fi J 

In the denominator, the quantity m/i — AqCos.^ is strictly positive as can be seen by 
either substituting Aq or comparing with the Newtonian value of (|36|) . By a power series 
expansion we get 



sin^ i-2fiE) 



1/2 



/i 



n i-2 



L'^{mijL — Aq cos^) 



^^.-2(0- 



(44) 



rr Aq 

Again, we can show that the numerator of the second term is proportional to sin^^. By 
using (p!3D, the second term of the function Hj{^) can be written as {mn — Aoy{mfi + AqY . 
For j > 0, we write Hj{^) as an algebraic sum of products of binomial expansions: 

^1 



1 

k,e=o 



-DH cos^+^ e + l±i-i£ cos^+2 e - 1 



(45) 



When k and i are of like parity, the square bracket in (^) has the form cos^'^ ^ — ^ hence 
we can take out a sin^ factor. When k is odd but i is even or conversely, the sum of pairs 
of terms {k, t) and (£, k) is proportional to cos^'''^"^ ^ — cos^'^'^'^^, with some integers qi and 
g2- This, again, contains the factor sin^^. Similarly, the corresponding expressions of Hj[^) 
for i = —1 and j = —2, given in Table 1, are proportional to sin^^. 

Thus we have shown that sin ^/rf is regular. 



(46) 



The integrand in (|39D , when written in terms of the the complex parameter 



w = e 



will contain a polynomial in w in the denominator. The possible roots lying inside the 
unit circle are at w = and at the root of rrifi — Aq cos ^ given in (|38| ) . Let us remark 
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that the factor m/i — AqCOS^ appears in the perturbative terms of Eq. (^) on the power 
n' + 2 — i = l — n — i, which is negative only for i > 2 — n. For these values of i the residue 
at wi is present. 

The function u in (|l]) this time can be a fraction with the numerator being a polynomial 
of sin C, and cos^ and the denominator an arbitrary power q > of mfi — AqCOsC,- In 



particular it may have the form (|33|) . The expressions for the azimuthal angle ip and the 



rate of change of the radius r, to be inserted in the perturbative term, are given by: 

mficos^-Ao . ( -2EV''^ Aqsiyi^ 

cos [iIj -ipo) = 2 ^= ~A c- ^^'^ 

m/i — COS 4 \ A* / mjj — AoCOSt, 

In the presence of uj the residue at wi appears if any of the conditions i > 2 — n , q > —n 
are satisfied. 

In summary, for the perturbed Keplerian motion (P), any integral @) with n < 
satisfying (P5D, by use of ( ^BD and ( ^7] ) can be turned into a complex integral in the 
variable w. Evaluation is simple since the residues TZ are in the origin and Wi {given in 



£ = [n{w = 0)+niw = wi)] . (48) 

If the radial equation (|]) is such that p < 2 — n and ip such that q < —n, then solely 
the first residue gives a contribution. In other cases, when the integrand f{w) has also a 
pole of order s at w = wi ^ oo its residue can be conveniently computed by: 

n(w = wi) = -^—T lim ^^[(w - WiYfiw)]. (49) 
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5. Computation of the period 



As a simple application of the method described in the previous section we compute 
the radial period for various perturbed Keplerian motions. 

The unperturbed system (yjj = 0) has the Kepler period 



\ i=0 / 

Here the corrections Tj from ipi may be obtained by putting n = — 2 in Eq. (p9D, substituting 
sin ^/rf from (0) and eliminating r by use of (0). Next, expressing r^^^ from (|16|), we get: 



By using the values of the functions A!j_ and Hi{C,) given in Table 1, we find that the 
correction Tj to the period vanishes for i = 2,3. In these cases, the poles of the period are 
in the origin of the complex w plane. For i = 4, the period is altered by a perturbative term 
given in the last column of the Table. 

Damour & Deruelle, 1985 compute the period T for perturbations of the form 
!fo = 3-E'^(3/i — m)/m(? and (fi = 2GE{7^ — 6m)/i/c^. (Damour and Deruelle enlist the 
constants A,...,H in an alphabetic notation; thus our energy E is distinct from their 
constant E). Their result can be recovered by summing up the corrections to for the 
entries i = and z = 1 in our Table 1 for c = 1. 




(50) 



The period T = Jq'^ of the perturbed system is the sum of contributions 




(51) 
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Table 1: Values of the constants A!|_, the function Hi{^) and corrections to the period T 
■ ^i~2 ^i~2 Hi_2{0 Correction Ti 

1 m/j2 /jAq m/i2 sin^ g yi 

EL2 _EL2 2_BL2 m/i-Ao cos £, 2m/t2 

2 2 -sin^^ 

3 2m/i 2y4o sin^ ^(Ag - m + "^Mo cos 

4 2(A2+mV) 4m/xAo sin^ e[2mV^o cos^ -|^(mV^-^o)^^^ 

-A2(A2+mV)cos2e-K-mV)'] 

6. Concluding Remarks 

We have analyzed in detail two types of generalized Keplerian parametrizations for 
which the evaluation of various integrals turns out to be a simple application of the residue 
theorem, Eqs. (|35| ) and ( ^) encapsulating the main results of this paper. We can state the 
generic result, that for a perturbed Keplerian motion (^) characterized by some value p > 0, 
the integrals (|l]) with n < min(l — p, — 1 — g) and n > 0, when parametrized by the eccentric 
and true anomaly parameters, respectively, are given by the corresponding residues at the 
origin TZ{0). For the evaluation of the remaining integrals with min(l — p, —1 — q) < n < , 
(such integrals can appear in the p > 3 case irrespective of the value of q), the computation 
of a second residue is necessary. We have done this in the eccentric parametrization, but a 
similar result holds for the true anomaly, where the second pole Zi = —Wi is the root of the 
equation + Aq cos x = 0. 

The solution of the radial motion we have presented has certain advantages over 
previous approaches ( [Rieth fc Schafer 1997] , |Ryan 19"96D and in our belief is at least a 



viable alternative to the use of the Laplace second integrals for the Legendre polynomials 
( IWhittaker fc Watson, 1996| , [Gopakumar fc Iyer, 199"7D . Although it is not a complete 



-17- 



parametrization of the motion [a review and comparison of such parametrizations ( |Brumbe"r^ 
19911 Pamour fc Deruelle, 1985| , [Blandford fc 'leukolsky, 1976HDamour fc Deruelle, 198q ) 



was given by |Khoner,S.A. fc Kopeikin,S.M., 1994|| , it is well suited when purely radial 



effects are considered. 

The power of the method we have described here was already tested in the evaluation 
of integrands involved in the radiative loss computations in [Gergely et al., 1998a| - |G"ergely| 
|et al., 1998c| . All integrands encountered there belong to the class n > 0, parametrizable 
by the true anomaly parameter, with the only pole in the origin. In this paper we have 
applied the eccentric anomaly parametrization for computing the period pertinent to a wide 
class of perturbations. A property holding for the spinning binary system, namely T = Tk 
was shown to be a consequence of the particular type of perturbing potential governing the 
motion, rather than a generic feature. 
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